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1 Introduction 

Because the running coupling g(T) of QCD decreases as a function of tem- 
perature T, the SU(3) color gauge theory becomes weakly coupled in the high 
temperature limit. Numerical simulations of the lattice gauge theory with dy- 
namical quarks have shown that the equation of state e(T) makes a transition 
from the low-temperature phase, where the effective degrees of freedom are 
hadrons, to a high temperature phase, where quarks and gluons can be excited 
independently. The transition temperature is T c w 150 MeV. 

Already in perturbation theory one finds that the quasiparticle excitations 
in the high-temperature phase are not massless. At low momenta, k < gT <C T 
(we are here considering the weak coupling limit), there are three propagating 
collective modes with frequency 

-W 2 k -^ul = \{l+N f /2)g*T\ 

the longitudinally and transversely polarized plasmons. At high momentum, 
the longitudinal mode is exponentially suppressed by its form factor; the trans- 
versely polarized modes carry an effective mass m* 2 = §iUp. 

This dynamically generated mass (self energy) is essential for the avoidance 
of infrared divergences which first occur at order g 4 in the equation of state. 
Incorporating the leading high-temperature contribution to the one gluon-loop 
diagram into effective gauge propagators and vertices, one obtains a consistent, 
i.e. gauge invariant, perturbative expansion of the partition function up to 
order g 5 \ng. Due to remaining infrared divergences attributable to static 
magnetic gauge interactions, the next term of order g 6 cannot be calculated 
perturbatively. 

Within the framework of perturbation theory, this divergence is seen to 
arise from the lack of a screening in the sector of time-independent transverse 
gauge fields at finite temperature. In a more formal context, the divergence 
is a result of the dimensional reduction of the high-temperature gauge theory, 
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which effectively becomes a three-dimensional theory as far as static quantities 
are concerned. The three-dimensional gauge theory is strongly coupled and 
confining, as revealed by the continued presence of an area law for space-like 
Wilson loops in the four-dimensional gauge theory at high temperature. In 
other words, the physics at length scales of order (g 2 T)~ x and larger remains 
nonperturbative and requires numerical calculations. 

2 Classical Color Plasma Description of Hard Thermal Loops 

The effective action describing the long-distance dynamics of the thermal gauge 
theory is well knownEfO'u 

r = r + Thtl = —ixF^F^ - 3u£ J ^teF^-HJ^f* (1) 

where = (1, v) is a light-like unit vector and denotes the gauge-covariant 
derivative. The second part of the action, Thtl, describes the contribution 
due to "hard thermal loops" , i.e. one-loop diagrams where the gauge bosons 
running around the loop carry at least thermal momenta (k > T). The action 
(1) generates effective n-point interactions with n > 2, in addition to the 
fundamental terms with n = 2, 3, 4. 

In spite of its formal simplicity, Thtl poses a number of non-trivial calcu- 
lational problems: 

• The effective action is highly nonlocal in space and time. 

• In Minkowski space the effective action acquires an imaginary 
part describing the damping of space-like fields. 

• To avoid double counting, soft momenta should be excluded from 
the hard thermal loops, but this should be done in a gauge invariant 
way. 

All these problems can be solved simultaneously by replacing Thtl with a ther- 
mal bath of classical, non-abelian charged particles that interact with the soft 
gauge fields. This is possible because the physics of hard thermal loops is ef- 
fectively classical: the thermal parts of loop diagrams are really tree diagrams, 
because the propagator in at least one section of the loop is on-shell. 

The equations describing the motion of classical, colored point particles in 
an external gauge field were first given by Wongul 

mi' 1 = p^ (2) 
mp* = gQ a F a ^p v (3) 
mQ a = -gf abc p»AlQ c . (4) 
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The second equation describes the effect of the non-Abelian Lorentz force, 
while the last equation describes the precession of the color charge vector as 
the particle moves through the field. This equation is really nothing else than 
an expression of the different definition of the color coordinate system at dif- 
ferent points in space. The equations do not describe effects of spin; these 
interactions are irrelevant in the HTL limit. If, instead of considering individ- 
ual point particles, one considers the evolution of a continuous distribution of 
such particles in phasa space, one obtains the non-Abelian transport equation 
first derived by Heinz El 



P" 



f(x,p,Q)=0, (5) 



which must be considered together with the Yang-Mills equation 

V" = ± / dpdQp v f(x >P ,Q)=r(x). (6) 



m t 

The particles are on-shell, i.e. f(x,p,Q) ~ 5{p 2 - m 2 )S(Q a Q a - Q 2 ). The 
transport equation (5) is gauge covariant, i.e. it remains unchanged under a 
gauge transformation 

UA^U- 1 - -Ud^U-\ (7) 
9 

f(x,p,Q)^ f(x,p,UQU- 1 ). (8) 

By considering the response of (5, 6) to a small external field A^, one can 
show that perturbations around a thermal distribution 

f eci (x,p,Q) = (e' 3 ^ -l)- 1 (9) 

lead to a set of equations that exactly mirror the HTL effective action (1). 
This can easily be understood in terms of the graphical representation of the 
one-loop contribution to the action of an ensemble of particles propagating in 
an external gauge field: 



Using p^ = mp^ and f(x,p, Q) = f e q(p) + Sf(x,p, Q), one finds 

D ab F b„u =g J dpdQ Q a v , 6f{xpj Q)= g J dp V »6f a (x,p) (10) 
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together with the explicit solution of (5) in the form 



oo 



Sf a (x,p) = -g 



d\ P \ J, 



duU ab (x, x-uv)- E b (x - uv), 



(11) 



where U ab (x,y) denotes the path-ordered link operator between x and y. The 
non-local nature of this solution is apparent. Retaining the differential form 
of (5) and solving this equation by means of test particles 



is more practical. In the limit of an infinite number of such test charges, 
obeying Wang's equations (2-4), (12) provides an exact solution of (5). 

One additional modification must be made in practice, in order to obtain 
a gauge covariant system of equations: the Yang-Mills equation (6) has to be 
formulated on a spatial lattice Q For the gauge fields, this technique is well 
known since the work of Wilson and Kogut and Susskind, the lattice definition 
of the color source term involves some subtleties. Essentially, the procedure is 
as follows: one surrounds each lattice site by a "cell" of sizes a 3 , where a is the 
lattice spacing, and associates the charge of every test particle momentarily 
contained in the cell with this lattice site. Currents are induced when test 
particles transit from one cell to the next; these lead to discontinuous charges 
of the electrical fields on the link orthogonal to the penetrated cell face. 

There are many technical details associated with the particle dynamics, 
which cannot all be discussed heretl Suffice it to say that electric fields act on 
the test particles only when they switch between cells, while magnetic fields 
act at each time step. If a particle does not have enough kinetic energy to 
transit into the next cell, it is reflected at the cell boundary. Careful ordering 
of the updates of the various quantities at each time step insures parity and 
time-reversal invariance, and an appropriate leap-frog algorithm conserves the 
phase-space measure. 

3 Results 

We have performed calculations on lattices as large as 32 3 with several hundred 
test particles per celllEm Energy conservation can be achieved at the level of 
10 -5 , and Gauss' law is conserved up to 10~ 10 . The parameters that can be 
freely chosen are: the gauge coupling g, the magnitude of the test particle color 
charge Q, the temperature T of their momenum distribution, and the average 
number of test particles per cell {no?). The HTL effective action contains these 




(12) 
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parameters only in the combination of the plasma frequency 



g 2 Q 2 n/T. 



(13) 



The independent variability of these parameters in the test particle formulation 
allows, in practice, for a numerical test of the correct scaling behavior. Let us 
now discuss some results. 
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Figure 1: Dispersion relation: free field 
in the continuum (dashed line: u) = k)\ 
free field on the lattice (dotted line: u>a = 
2 sin (ka/2)); plasma modes: transverse 
(theory: the upper solid curve; numerical 
results: solid rectangles) and longitudinal 
(theory: the lower solid curve; numerical 
results: open rectangles). [From reflj] 



The dispersion relation obtained on a 16 3 lattice is shown in Fig. 1. At 
low wave numbers k, the eigenfrequencies nicely follow the expectation from 
continuum HTL perturbation theory (solid lines), whereas at large fc, the fre- 
quencies approach the free lattice dispersion curve (dotted line). For k ^ 0, 
longitudinal and transverse modes clearly have different frequencies, and the 
plasma frequency w p is obtained for k = 0. The dependence of w p on the 
combination gQ\JnjT agrees very nicely with the HTL prediction, as shown 
in Figure 2. The non-Abelian plasma modes are strongly damped. Again, 
the analysis shows that the plasmon damping rate 7 P for k — is in excellent 
agreement with HTL perturbation theory (see Fig. 3). 

While this is encouraging, it is important for many applications to verify 
the correct HTL behavior of the response functon for arbitrary combinations 
of lo and k. In the Abelian limit, this can be achieved by a linear response 
analysis, when the system of gauge field and test particles is driven by a weak 
current of the form 



t) = j jo sin(urf) sin (k ■ x \ . 



(14) 



The dissipative and dispersive part of the response function can then be ob- 
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Figure 2: Space-time average of j ■ E plot- 
ted as a function of uj/k (longitudinal polar- 
ization). Upper window: in-phase average; 
lower window: out-of-phase average. [From 



Figure 3: Space-time average of j ■ E plot- 
ted as a function of uj/k (transverse polar- 
ization). Upper window: in-phase average; 
lower window: out-of-phase average. [From 
refT 



tained by means of the in-phase and out-of-phase averages: 

(j(x,t) ■ E(x,t)) and (j(x, t) ■ E(x, t - ^)), (15) 

respectively. Transverse and longitudinal response are simply obtained by 
choosing j _L k or j || k. Again, one finds nice agreement with the HTL 
predictions for both response functions, as shown in Figs. 4, 5. The spikes 
of the response functions at the plasma resonance frequency cu p (k) are clearly 
visible. The smallness of the dispersive part of the transverse response for 
u> < k is indicative of the absence of magnetic screening at the HTL level. 

As an application, we have studied the diffusion rate of the Chern-Simons 
number A^cs m the unbroken 5(7(2)-Higgs Theory, which corresponds to the 
high-temperature phase of the electroweak theory in the limit sin 2 9\y — * 0. 
It was long suspected, on the basis of scaling arguments that this rate Tcs 
should be proportional to (g 2 T) A . However, Arnold, Son, and Yaffe argued 
convincingly that the rate should be suppressed by another factor g 2 , because 
space-like modes of the gauge field are strongly damped by the HTL dynamics. 
Since 

H H TL(fc, W ) ^ = (16) 

the relevant frequency domain for long-distance fluctuations of the gauge field 
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Figure 4: Scaling behavior of the plasma 
frequency. Solid line: lu p = 0.39gQy^ {n)/T 
(theory); trianglps: numerically obtained 
values. [From refB] 



Figure 5: Scaling behavior of the plasma 
rates. Solid line: 7 P = 0.176g 2 T (the- 
ory); triangles* numerically obtained val- 
ues. [From ref|3| 



becomes 

u~— 2 g^T (for k~g 2 T). (17) 

Our numerical simulations of Tcs show that, indeed, the rate depends linearly 
on l/u)p as expected from (17), confirming the ASY scaling law (see Fig. 6)]j 



4 Summary 

The infrared limit of thermal gauge field dynamics in Minkowski space can be 
represented by the local interaction of the classical gauge field with a thermal 
plasma of classical colored particles. A lattice formulation of this dynamical 
system can be derived which conserves all important symmetries and lends 
itself to numerical simulation. We have verified that the HTL dispersion re- 
lations known from perturbation theory are well reproduced, and we have 
confirmed the Arnold-Son- Yaffe scaling relation for the Chern-Simons number 
diffusion rate in the Yang-Mills-Higgs model. 

Two interesting generalizations for future work are: The study of the evo- 
lution of gauge field configurations far off equilibrium, such as they occur in 

a More recently, Bodecker has argued that the dissipative response of the hard thermal loops 
has logarithmic corrections ~ log(l/g). We observed no indication for such a behavior, but 
the logarithm may be difficult to see in the numerical results. 
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Figure 6: Result for k at three values of 
average particle number (n) but a common 
value of lattice spacing and other physi- 
cal parameters (g, Q, T), plotted against 
9 4 T 2 /m^. The error bars in mi. reflect 
uncertainty in the damping from hard clas- 
sical lattice modes. The old picture predicts 
a flat line, while the ASY picture predicts 
a straight line through Wie origin, like the 
illustrated fit. [From refl 



relativistic heavy ion reactions, and the inclusion of hard interactions among 
the colored particles, generalizing the non-Abelian Vlasov equation to a non- 
Abelian Boltzmann equation and effectively introducing a color mean field into 
the parton cascade modelt2l 
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